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1 Introduction 

In this paper, we consider a class of nonlocal feedback controllers with appli¬ 
cation to the control of certain nonlinear partial differential equations. The 
research on feedback control laws of this type has become quite active in theo¬ 
retical physics for stabilizing wave-type solutions of reaction-diffusion systems 
such as the Schlogl model (also known as Nagumo or Chafee-Infante equation) 
or the FitzHugh-Nagumo system. 

The controllers can be characterized as follows: First of all, they are a 
generalization of Pyragas type controllers that became very popular in the 
past. We refer to 120], |2l], and the survey volume [24]. In the simplest form of 
Pyragas type feedback control, the difference of the current state u(a:, t) and 
the retarded state u{x,t — r), multiplied with a real number k, is taken as 
control, i.e. the feedback control / is 


f{x, t) := K (u(x, t) — u{x, t — t)). 


where r is a fixed time delay and k is the feedback gain. 

In the nonlocal generalization we consider in this paper, the feedback con¬ 
trol is set up by an integral operator of the form 



( 1 ) 


Here, different time delays appear in a distributed way. Depending on the par¬ 
ticular choice of the kernel g, various spatio-temporal patterns of the controlled 
solution u can be achieved. We refer to Ejdaiis], and |25| with application to 
the Schlogl model and to niiiniiss] with respect to control of ordinary differ¬ 
ential equations. 

Our main goal is the selection of the kernel g in an optimal way. We want to 
achieve a desired spatio-temporal pattern for the resulting state function and 
look for an optimal feedback kernel g to approximate this pattern as closely as 
possible. For this purpose, in the second half of the paper we will concentrate 
on a particular choice of g as a step function. 

We are optimizing feedback controllers but we shall apply methods of op¬ 
timal control to achieve our goal. This leads to new optimal control problems 
for reaction-diffusion equations containing nonlocal terms with time delay in 
the state equation. We develop the associated necessary optimality conditions 
and discuss numerical approaches for solving the problems posed. Working 
on this class of problems, we observed that standard quadratic tracking type 
objective functionals are possibly not the right tool for approximating desired 
time-periodic patterns. We found out that the so-called cross correlation par¬ 
tially better fits to our goals. We report on our numerical tests at the end of 
this paper. 

This research contributes results to the optimal control of nonlinear re¬ 
action diffusion equations, where wave type solutions such as traveling wave 
fronts or spiral waves occur in unbounded domains. We mention the papers 
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mini on the optimal control of systems that develop spiral waves or dsuni 
115] on systems with heart medicine as background. Moreover, we refer to j5l 
[7] , where different numerical and theoretical aspects of optimal control of the 
Schlogl or FitzHugh-Nagumo equations are discussed. It is a characteristic fea¬ 
ture of such systems that the computed optimal solutions might be unstable 
with respect to perturbations in the data, in particular initial data. 

Feedback control aims at generating stable solutions. Various techniques 
of feedback control are known, we refer only to the monographies [lOlITTlfTSl 

and to the references cited therein. Moreover, we mention m on feedback 
stabilization for the Schlogl model. Pyragas type feedback control is one as¬ 
sociated field of research that became very active, cf. |M] for an account on 
current research in this field. In associated publications, the feedback control 
laws were considered as given. For instance, the kernel in nonlocal delayed feed¬ 
back was given and it was studied what kind of patterns arise from different 
choices of the kernel. 

The novelty of our paper is that we study an associated inverse (say de¬ 
sign) problem: Find a kernel such that the associated feedback solution best 
approximates a desired pattern. 


2 Two models of feedback control 


We consider the following semilinear parabolic equation with reaction term R 
and control function (forcing) /, 

dtu — Au + R{u) = f (I) 

subject to appropriate initial and boundary conditions in a spatio-temporal 
domain 17 x {0,T). Using a feedback control in the form 0 , we arrive at the 
following nonlinear initial-boundary value problem that includes a nonlocal 
term with time delay, 

g{T)u{x, t — t) dr — u{x, t)^ in 17, 

u(x, s) = uo(x, s) in 17, 

dnu{x, t) = 0 on r, 

( 2 ) 

for almost all t G (0,T), s G [—T, 0]. 

Here, denotes the outward normal derivative on F = dfl. We want to 
determine a feedback kernel g G L°°{Q,T) such that the solution u to 0 is 
as close as possible to a desired function Ud- The function g will have to obey 
certain restrictions, namely 


dtu(x, t) — Au{x, t) + R{u{x, t)) = K 


0 < git) < p 
/ gis)ds = l, 


a.e. on [0,r], 


0 


(3) 

(4) 
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where /? > 0 is a given (large) positive constant. This upper bound is chosen 
to have a uniform bound for g. It is needed for proving the solvability of the 
optimal control problem. 

We shall present the main part of our theory for the general type of g 
defined above. In our numerical computations, however, we will concentrate 
on functions g of the following particular form: We select ti, t 2 such that 
Q <ti < t 2 <T, t 2 — ti > 5 > Q and define 


9{i) 


1 

-, ti<t<t2 

t2-ii " 

0, elsewhere. 


(5) 


It is obvious that g satisfies the constraints @.0 with 13 = 1/5. Using this 
form for g, we end up with the particular feedback equation 


dtu{x, t) — Au{x, t) + R{u{x, t)) = K (- / u{x,t — t) dr — u{x, t) ) . 

V^2 - h Jti J 

( 6 ) 


In ^ , we will also vary k in the state equation as part of the control variables 
to be optimized. In contrast to this, k is assumed to be fixed in the model with a 
general control function g. In the special model, we have a restricted flexibility 
in the optimization, because only the real numbers k, ti, t 2 can be varied. Yet, 
we are able to generate a class of interesting time-periodic patterns. 

Throughout the paper we will rely on the following 


Assumptions. The set f2 C TV < 3, is a bounded Lipschitz domain; 
for A = 1, we set 17 = (a, 5). By T > 0, a finite terminal time is fixed. In 
theoretical physics, also the choice T = oo is of interest. However, we do not 
investigate the associated analysis, because an infinite time interval requires 
the use of more complicated function spaces. Moreover, the restriction to a 
bounded interval fits better to the numerical computations. Throughout the 
paper, we use the notation Q := il x (0,T) and S = F x (0,T). for the 
space-time cylinder. 


Remark 1 We will often use the term ’’wave type solution” or ’’traveling wave”. 
This is a function (x, t) i—>■ u{x, t) that can be represented in the form u{x, t) = 
v{x — ct) with some other smooth function v. Here, c is the velocity of the 
wave type solution. Such solutions are known to exist in 17 = K but not in in 
a bounded interval 17 = (a, h). 

In our paper, the terms ” wave type solution” or ’’traveling wave” stand 
for solutions of the Schlogl model in the bounded domain (a, h). We use these 
terms, since the computed solutions exhibit a similar behavior as associated 
solutions in 17 = K. 


The reaction term R is defined by 

R{u) = p{u-ui){u-U 2 ){u-U 3 ,), (7) 

where ui < U 2 < and p > 0 are fixed real numbers. In our computational 
examples, we will take p := 1. The numbers f = 1,..., 3, define the fixed 
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points of the (uncontrolled) Schlogl model Q. In view of the time delay, we 
have to provide initial values uq for u in the interval [—T, 0] for the general 
model and in [—^ 2 ) 0 ] for the special model We assume uq g C(f2 x 
[—T, 0]) or uo G C(f2 x [—^2,0]), respectively. The desired state Ud is assumed 
to be bounded and measurable on Q. 


3 Well-posedness of the feedback equation 

In this section, we prove the existence and uniqueness of a solution to the 
general feedback equation To this aim, we first reduce the equation to an 
inhomogeneous initial-boundary value problem. For t G [0,7’], we write 

/.T .t .T 

/ g{T)u{x,t — t) dr = / g{T)u{x,t — t) dr + / g{T)u{x,t — t) dr 

Jo Jo Jt 


= -.Ug{x,t) 


= / g{T)u{x,t — t) dr + Ug{x,t). 

Jo 


The function Ug is associated with the fixed initial function uq and is defined 
by 


Ug{x,t) = j^ g{T)uo{x,t-T) dr- 


notice that we have t — t < 0 in the integral above. By the assumed continuity 
of Uo, the function Ug belongs to C(l7 x [0, Tj). 

Next, for given g G T^(0, T), we introduce a linear integral operator K{g) : 
L^iQ) -G L\Q) by 

{K{g)u){x,t)-.= [ g{T)u{x,t- T)dT. ( 8 ) 

Jo 

Substituting s = t — r, -we obtain the equivalent representation 

{K{g)u){x,t) = / g{t- s)u{x,s)ds. 

^0 

Inserting Ug and K{g) in the state equation we obtain the following non¬ 
local initial-boundary value problem: 


dtu — Au + R{u) + Ku — K K{g)u = nUg in Q, 

u(x,0) = uo(x,0) in 12, 
dnU = 0 on T". 

In the next theorem, we use the Sobolev space 

W(0,T) = L^(0,T;II^(f2)) n Id^(0,T; L^(f2)). 


( 9 ) 
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Theorem 1 For all g € L°°{0,T), Ug € L^(Q), p > |, and uq € C(i7 x 
[—T, 0 ]), the problem has a unique solution u G TT(0,T) r\C{Q). 

Proof We use the same technique that was applied in | 8 ] to show the existence 
and continuity of the solution to the FitzHugh-Nagumo system. Let us mention 
the main steps. First, we apply a simple transformation that is well-known in 
the theory of evolution equations. We set 

u = 


with some A > 0. This transforms the partial differential equation in © to an 
equation for the new unknown function v, 

Vt — Av F -I- (A -I- k)v = KKx{g)v + e~^*KUg, (10) 


where the integral operator K\{g) is defined by 

{Kx{g)v){x,t) = f s)v{x,s)ds. 

Jo 

li g G L°°{0,T), then both operators K{g) and K\(g) are continuous linear 
operators in Lp{Q), for all p > 1. Moreover, due to the factor the 

norm of K\{g) : Lf{Q) -G L^{Q) tends to zero as A —?> oo. We obtain 

\\K\{g)\\c(L^(Q)) < ;^I|5||l“(o.t) (H) 

with some constant c > 0. To have this estimate, we assumed in ^ that g is 
uniformly bounded by the constant /3. If A is sufficiently large, then we have 

f [e~^*R{e^*'v) + {X +k)v - KKxig)v]vdxdt > ^\\v\\'l2tQ\ Vu G L^(Q), 

JQ ^ 


because the coercive term {X + k)v in the left side is dominating the other 
terms, cf. j^. 

With this inequality, an a priori estimate can be derived in Lf{Q) for any 
solution V of the equation ([^. Now, we can proceed as in [ 8 ]: A fixed-point 
principle is applied in Lf{Q) to prove the existence and uniqueness of the 
solution V that in turn implies the same for u. For the details, the reader 
is referred to |S], proof of Theorem 2.1. However, we mention one important 
idea: Thanks to (111, the term (A-|-k) absorbes the non-monotone terms in the 
equation (101 so that, in estimations, equation ( 10 ) behaves like the parabolic 
equation 

Vt — Av + R{v) = F 


with a monotone non-decreasing nonlinearity R and given right-hand side 
F G L^(Q), p > 5/2. This fact can be exploited to verify, for each r > 0, 
the existence of a constant C^. > 0 with the following property: If 5 G L°°{Q) 
obeys ||g||L~(Q) < r and u is the associated solution to (|^, then 


||u||l~(Q) < Cr- 


( 12 ) 


□ 
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4 Analysis of optimization problems for feedback controllers 

4.1 Definition of two optimization problems 
General kernel as eontrol 


Let a desired function Ud € L°°{Q) be given. In our later applications, Ud 
models a desired spatio-temporal pattern. Moreover, we fix a non-negative 
function cq G L°° (Q) . This function is used for selecting a desired observation 
domain. We consider the feedback equation ([^ and want to find a kernel g 
such that the associated solution u approximates Ud as close as possible in 
the domain of observation. This goal is expressed by the following functional 
j : L^{Q) X L°°{0,T) -G M that is to be minimized. 





V 

2 





Here, zz > 0 is a Tikhonov regularization parameter. The standard choice of cq 
is CQ{x,t) = 1 for all {x,t) G Q. Another selection will be applied for periodic 
functions Ud'- c{x, t) = 1 for all {x, t) G Q with t>T/2 and c(a;, t) = 0 for all 
{x,t) G Q with t <Tf2. 

By Theorem]^ to each g G L°°(0, T) there exists a unique associated state 
function u that will be denoted by Ug. Then j does only depend on g and we 
obtain the reduced objective functional J, 


J -9^ j{ug,g)- 

Therefore, our general optimization problem can be formulated as follows: 


mm. 
sec 


Ji 9 )-=^JJ^CQ{ug-Udfdxdt+^J g'^{t)dt, (PG) 


where C C L°°{Q,T) is the convex and closed set defined by 

C := € L°°{0,T) : 0 < g{t) < /3 a.e. in [0,T] and j g(t) dt = l.| 


Notice that C is a weakly compact subset of L^(0,T). The restrictions on g 
are motivated by the background in mathematical physics. In particular, the 
restriction on the integral of g guarantees that 

g(T)u{x, t — T)dT — u(x, t) = 0, 

if u{x,t — T) = u{x, t) in Q. By the definition of Ug , the optimization is subject 
to the state equation (§. 
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Special kernel as control 

The other optimization problem we are interested in, uses the particular form 
^ of the kernel g, 


min Js(K,ti,t2) 



where is the solution of for a given triplet (k, ti, ^ 2 )- This problem 

might fail to have an optimal solution, because the set of admissible triplets 
{K,ti,t 2 ) is not closed. Notice that we need ti < t 2 in Therefore, we fix 
i5 > 0 and define the slightly changed admissible set 


Cs := {(k, ti,t2) G : 0 < ti < t2 < T, t2 — ti > S, k G M} 


that is compact. In this way, we obtain the special finite-dimensional optimiza¬ 
tion problem for step functions g, 

min Js{K,ti,t 2 ) ■■= 1: CQ{u^n^ti,t 2 )-Ud)'^dxdt+’^{tl+tl +K^). 

^JJo ^ 

(PS) 


4.2 Discussion of (PG) 

The control-to-state mapping G 

Next, we discuss the differentiability of the control-to-state mappings g 1 -^ Ug 
and {K,ti,t 2 ) I—>■ u{K,ti,t 2 ). First, we consider the case of the general kernel 
g. The analysis for the particular kernel (|^ is fairly analogous but cannot 
deduced as a particular case of (PG). We will briefly sketch it in a separate 
section. 

By Theorem we know that the mapping G : g ^ Ug is well defined 
from L°°(0,T) to C{Q). Now we discuss the differentiability of G. To slightly 
simplify the notation, we introduce an operator K : L°°{0,T) x G{Q) —>■ G{Q) 
by 


lC{g,u) = K{g)u, 

where K{g) was introduced in ([^; notice that /C is bilinear. Let us first show 
the differentiability for K,. 
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We fix g G L°°(0, T), u G C{Q), and select varying increments h G L°°(0, T), 
V G C(Q). Then we have 

IC{g + h,u + v)= ( [g{T) + h{T)][u{x,t — t) -\-v{x,t — T)]dT 

Jo 

= / g{T)u{x,t — t) dr + / h{T)u{x,t — t) dr + / g{T)v{x,t — t) dr 
Jo Jo Jo 

' -V-^ 

A{g,u)(h,v) 

+ / h{T)v{x,t — t) dr = lC{g,u) + A{g,u){h,v) + R{h,v), 

Jo 

■' -V-" 

R(h,v) 


where A{g,u) : L°°{0,T) x C{Q) -G C{Q) is a linear continuous operator and 
R : L°°{0,T) X C(Q) -G C{Q) is a remainder term. It is easy to confirm that 


\\R{h, ^')llc(Q) 

||(/l,ti)||j;,oo(o ^)xC(Q) 


if ll(^!'*^)llL“(o,r)xC(Q) 0- 


Therefore, JC is Frechet-differentiable. As a continuous bilinear form, JC is also 
of class 

Now, we investigate the control-to-state mapping G : L°°{0,T) -G C{Q) 
defined hy G : g ^ Ug, where the state function Ug is defined as the unique 
solution to 


dtu — Au + R{u) + Ku = KIC{g,u) + nUg in Q 

dnU = 0 in H (13) 

u(0) = Uo(0) in G. 

In what follows, the initial function uq will be kept fixed and is therefore not 
mentioned. Of course, Ug, G and some of the operators below depend on uq, 
but we will not explicitely mention this dependence. To discuss G, we need 
known properties of the following auxiliary mapping Q : v i-A u, where 

dtU — Au + R{u) + KU = V in Q 
dnU = 0 in A 
m( 0) = Mo(0) in J7. 


This mapping Q is of class from LP{Q) to W{0,T) n C(Q), if p > |, in 
particular from L°°{Q) to L°°{Q), cf. |S] or, for monotone R, |^, |53], |28| . 

Now (consider v := k {lC{g,u) + Ug) as given and keep the initial function 
Uq fixed), u solves (13l if and only if u = Q{KlC{g,u) + KUg), i.e. 


u — Q{KK,{g,u) + uUg) = 0. 

We introduce a new mapping F : L°°{Q) x L°°{Q,T) -G L° 


( 14 ) 

’ (Q) defined by 


F{u,g) := u — G{KlC{g,u) + uUg). 
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Then, (141 is equivalent to the equation 

T{u,g) = 0. 


(15) 


We have proved above that the mapping {g,u) JC{g,u) is of class from 
L°°(0,r) X L°°{Q) to L°°{Q). Obviously, also the linear mapping g Ug is 
of class from L°°(0,T) to L°°{Q). By the chain rule, also T is from 
L°°{Q) X L°°(0,T) — >• L°°{Q) and the mappings dgT{u,g), duJ^{u,g) are 
continuous in the associated pairs of spaces. 

To use the implicit function theorem, we prove that du^{u,g) is continu¬ 
ously invertible at any fixed pair {u,g). Therefore, we consider the equation 


duT{u,g)v = z 


(16) 


with given right-hand side z € L°°{Q) and show the existence of a unique 
solution V G L°°{Q). The equation is equivalent with 

V — Q'{KlC{g,u) + KUg)KlC{g,v) = z. (17) 

P 

Writing for convenience p = KlC{g,u) + kUq, we obtain the simpler form 

V — Q'{p)K, K{g)v = z. 


A function z € L°°{Q) does not in general belong to W(0,T). To overcome 
this difficulty, we set w := v — z and transform the equation to 

w = Q'{p) K,K{p){w + z) = Q'{p)q. (18) 

Q 

where q := KK{g){w + z). As the next result shows, w is the solution of a 
parabolic PDE, hence w £ W(0,T). 

Lemma 41 Let q £ L^^Q) with p > hl‘1 he given. Then we have y = Q'{p)q 
if and only if y solves 


dty - Ay + R'{u)y + Ky = q in Q 
dny = 0 in E 
2 /( 0 ) =0 in f2, 

where u is the solution associated with p, i.e. 

dtu — Au -I- R{u) + Ku = p in Q 
dnU = 0 in E 
u(0) = uq in Q. 
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We refer to [S]. For monotone non-decreasing functions R, this result is well 
known in the theory of semilinear parabolic control problems, see e.g. m,m, 
or [251 Thm. 5.9]. By Lemma 41 the solution w of (18 1 is the unique solution 
of the linear PDF 


{dtw — Aw -I- R'{u)w + Kw){x,t) = q{x,t) 


= K / g{T)w{x,t — t) dr + K / g{T)z{x,t — t) dr 

Jo Jo 


(19) 


subject to ■u;(0) = 0 and homogeneous Neumann boundary conditions. By the 
same methods as above we find that, for all z G L°°{Q), equation (19) has a 
unique solution w G W(0,T) n L°°{Q). 

After transforming back by u = w -I- z, we have found that for all z G 
L°°{Q), (16 1 has a unique solution v G L°°{Q) given hy v = w + z. Therefore, 
the inverse operator duJ^{u,g)~^ exists. The continuity of this inverse map¬ 
ping follows from a result of |5] that the mapping z i-G w defined by (191 is 
continuous in L°°{Q). 

Next, we consider the operator dgfF. It exists by the chain rule and admits 
the form 

dgfF{u,g)h = Q'{k {K{g)u + Ug))K {K{h)u + dgUgh). 

Setting again p = k {K{g)u Ug) and q = k {K{h)u + dgUgh), we see that 

dgfF{u,g)h = f], 

where, by Lemma rj solves the equation 

dtT] — Arj + R'{u)r] + ng = q = KK{h)u + ndgUgh 

subject to homogeneous initial and boundary conditions. Therefore, g is the 
unique solution to 

{dtg — Ag + R'{u)g + Kg){x,t) = K / h{T)u{x,t — t) dr 

Jo 

r 

+K / h{T)uo{x,t — t) dr 

g{x,0) = 0 
dng = 0 . 

By u{x, t) = uo(x, t) for —T < t < 0, we can re-write this as 

{dtg — Ag + R'{u)g + Kg){x,t) = K / h{T)u{x,t — t) dr 

Jo 

dng = 0 
g{x, 0 ) = 0 . 

Again, the mapping h i-G w is continuous from L°°{0,T) to hF(0,T) r)C{Q). 
Collecting the last results, we have the following theorem: 
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Theorem 2 (Differentiability of G) The control-to-state mapping G : g <—>■ 
Ug assoeiated with equation © is of elass G^. The first order derivative z := 
G'{g)h is obtained as the unique solution to 

r 

{dtz — Az + R'{ug)z + K z){x,t) = K / h{T)ug(x,t — t) dr 

Jo ^ 

/ g {T)z{x,t — t) dr in Q 

Jo 

dnZ = 0 in S 

z{-,t) = 0, —T <t<0 in f2. 

( 20 ) 

Proof We already know by Theorem]^ that, for all g G L°°{0,T), there exists 
a unique solution u = G{g) G Wfid,T) r\C{Q) solving the equation 

Riu,g) = 0 . 


We discussed above that the assumptions of the implicit function theorem are 
satisfied. Now this theorem yields that the mapping g i—)■ G{g) is of class C^. 

The derivative G'{g)h is obtained by implicit differentiation. By definition 
of G{g), we have 


{dtG{g) - AG{g) + R{G{g)) + KG{g))ixfi) = K [ gir) G{g)ixfi - t) dr 

Jo 

r 

+K / g{T) uo{x,t — t) dr 


dnG{g) = 0 

Gig)i-,t) = uoi-,t), -T <t<0. 


( 21 ) 


Implicit differentiation yields that z := G'{g)h is the unique solution of ( [20| . 
Notice that 

/ g{T)G{g){x,t-T)dT+ j g{T)uo{x,t-T)dT = g^r) G{g){x,t-T) dr. 
Jo Jt Jo 

□ 


4.3 Existence of an optimal kernel 

Theorem 3 For all n > 0, (PG) has at least one optimal solution g. 

Proof Let (gn) with gn G G for all n G N be a minimizing sequence. Since C is 
bounded, convex, and closed in L°°(0,T), we can assume without limitation 
of generality that gn converges weakly in L^{0,T) to g, i.e. gn ^ g, n ^ oo. 
The associated sequence of states obeys the equations 

dtUn - AUn + KUn = dn ■■= - K R{Un) + KK{gn)Un + nUg. ( 22 ) 
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By the principle of superposition, we split the functions u„ as = it + 
Un 7 where u is the solution of (221 with right-hand side (i„ := 0 and initial 
value it(0) = ito(O), while it„ is the solution to the right-hand side defined 
above and zero initial value. In view of (121, all state functions ii„, hence also 
the functions are uniformly bounded in L°°{Q). Thanks to |11[ Thm. 4], 
the sequence (i(„) is bounded in some Holder space By the Arzela- 

Ascoli theorem, we can assume (selecting a subsequence, if necessary) that Un 
converges strongly in L°°{Q). Adding to the fixed function ii, we have that 
(it„) converges strongly to some u in L°°{Q). 

The boundedness of (ii„) also induces the boundedness of the sequence ((i„) 
in L°°{Q), in particular in L^{Q). Therefore, we can assume that converges 
weakly in to d, n —>■ oo. Since (un) is the sequence of solutions to the 

"linear” equation ( p^ with right-hand side d„, the weak convergence of (dn) 
induces the weak convergence of ii„ ^ it in W{0,T), where u solves (22) with 
right-hand side d. 

Finally, we show that 

d{t) = —KR{u{t)) + K {K{g)u){t) K Ug{t) 

so that u is the state associated with g. Obviously, it suffices to prove that 
K{gn)un converges weakly to K(g)u in L^(Q). To this aim, let an arbitrary 
(fi G L‘^{Q) be given. Then we have 


JJ gn{T)un{x,t - t) dr^ dxdt 

= j 9n{'r){^j J^ip{x,t)Unix,t-T) 


— t) dtdx dr. 


(23) 


Clearly, the strong convergence of (it„) in L°°{Q) yields 


ip{x, t)Unix, t — ■) dtdx 


g}{x, t)u{x, t — ■) dtdx 


in L^{0,T). Along with the weak convergence of g^, this implies 


lim 

n—>-oo 


9u{t) 


ip{x, t)un{x, t — t) dtdxdr 


= / 9i^) / / ‘fix, t)u{x, t — t) dtdxdr. 

Jo J T J Q 


In view of (231, we finally arrive at 


// (p{x,t) / gn{T)un{x,t—T) drdxdt ^ // <g{x,t) / g{T)u{x,t—T) drdxdt 

JJq Jo JJq Jo 

as n —>■ oo. Since this holds for arbitrary (p G L^{Q), this is equivalent to the 
desired weak convergence K{gn)un K{g)u in L'^{Q). □ 
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4.4 Necessary optimality conditions 
4-4-^ Adjoint equation 

In the next step of our analysis, we establish the necessary optimality condi¬ 
tions for a (local) solution g of the optimization problem (PG). This optimiza¬ 
tion problem is defined by 


min J(g), 

0 < g(t) < j3 for almost all t G [0, T], 


[ g{T)dT = l. 

Jo 


(24) 


Although the admissible set belongs to L°°{0,T), we consider this as an opti¬ 
mization problem in the Hilbert space L^(0,T). 

To set up associated necessary optimality conditions for an optimal solution 
of (241, we first determine a useful expression for the derivative of the objective 
functional J. We have 

= CQ{ug - Udf dxdt + J g{tf dt 

= CQ{G{g) - Ud)'^ dxdt + J g{tf dt. 


Let now be an arbitrary (i.e. not necessarily optimal) g G L°°(0,T) be given 
and let u = G{g) be the associated state. Then we obtain for h G L°°(0,T) 

J'{g)h = v f g{t)h{t)dt+ ff cq{u — Ud){G' (v)h) dxdt 
Jo J JQ 

= / vg{t)h{t)dt-G // CQ{x,t){u{x,t) — Ud{x,t))z{x,t)dxdt (25) 

Jo J Jo 


with the solution z to the equation (201 for Ug := Ug = u. 

The implicit appearance of h via z can be converted to an explicit one by 
an adjoint equation. This is the following equation: 

r 

{—dtip — Glip + R'{u)if + K^){x,t) = K / g{T)(p{x,t + t) dr 

Jo 


+cq{x, t){u{x, t) - Ud{x, t)) 


a.e. in Q, 

dn^p = 0 in 27, 

V5(-,t)=0 te[T,2T]. 


(26) 


The solution tp of (261 is said to be the adjoint state associated with g. 
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Lemma 1 Let g, h € L°°(0,T), and u = Ug be given. If z is the solution to 
the linearized equation (201 for Ug := u and (p is the unique solution to the 
adjoint equation (26 I, then the identity 



u — Ud) z:)(x, t) dxdt = k 



h{T)u{x, t — t) dr 


is fulfilled: 


dxdt 

(27) 


Proof We multiply the first equation in (20 1 by the adjoint state (p as test 
function and the first equation in (26) by z. After integration on Q and some 
partial integration with respect to x, we obtain 


JJ {dtz p + Vz ■ Vp + {R'{u) + k)z p) dxdt 

= K JJ g{T)z{x,t — t) dr^ p{x,t) dxdt 
+ K JJ ij h{T)u{x,t — T)dTjp{x,t)dxdt 


and 


JJ {—z dtp + Vz ■ Vp + {R'(u) + k)z p) dxdt 
= K JJ iJ g{T)p{x,t + T)dTjz{x,t)dxdt 


Integrating by parts with respect to t, we see that 


cq{u — Ud) z dxdt. 


{—z)dtpdxdt= // pdtzdxdt] 

Q J Jq 


notice that we have z(0) = 0 and p{T) = 0. Comparing both weak formulations 
above, it turns out that we only have to confirm the equation 


// / g{T)z(x,t—T) p{x,t) dr dxdt = // / g{T)p{x,t+T) z{x,t) dr dxdt. 

JJqJo JJqJo 

( 28 ) 
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Then the claim of the Lemma follows. To show (281 , we proceed as follows: 
rT 

IQ Jo 


// / g{T)z{x,t — T)(p{x,t) dr dxdt 

J Jq Jo 

= / g{T) z{x,t — t) (p{x,t) drdtdx 

Jn Jo Jo 

= / / / g{T) z{x,t — t) (p{x,t) drdtdx 

Jq Jo Jo 

= / / / g{t — 7]) z{x,ri) <f{x,t) dgdt dx 

Jq Jo Jo 

= / g{t — rf) (p{x,t) dt z{x,g) dg dx 

J J 0 J rj 

= / g{a) (p{x,g + a) da z{x,g) dg dx 

Jq Jo Jo 

= / g{a) (p{x,g + a) da z{x,g) dg dx 

Jq Jo Jo 


(29) 


g{T) (p{x, t + t) z{x, t) dr dxdt. 


Q Jo 


We used z{x,t — r) = 0 for r > t in the second equation, the substitution 
g = t — T in the third, the Fubini theorem in the fourth, the substitution 
a = t — g in the fifth, the property <f{x, t) = 0 for t >T in the sixth equation. 
Finally, we re-named the variables. □ 

Corollary 1 At any g G L°°{0,T), the derivative J'{g)h in the direetion 
h G L°°{0,T) is given by 

J'{g)h = J ly g{t) h{t) dt + K J (^JJ ^ ~ 

where Cp is the unique solution of the adjoint equation (26l. 


This follows immediately by inserting the right-hand side of (27) in (25) 
and by interchanging the order of integration with respect to t and t. 


4..4..2 Neeessary optimality eonditions for (PG) 

Let us now establish the necessary optimality conditions for an optimal solu¬ 


tion g ol (24). They can be derived by the Lagrangian function L : L°°(0, T) x 


( [ - 1 ) • 


L{g,fi) := J{g)+fJ. 


0 
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If g is an optimal solution, then there exists a real Lagrange multiplier fi such 
that the variational inequality 

J'i9){9 -9) + P‘[ { 9 - 9 )dt >0 for all 5 > 0 
Jo 

is satisfied. Inserting the result of Corollary [l] for h := g — g, we find 

J ^i^g{t) + p. + K jj (p{x, s)u{x, s — t) dxds^ {g{t) — g{t)) dt > 0 (30) 

for all 0 < g < /3. 

Remark 2 For a Lagrange multiplier rule to hold, a regularity condition must 
be fulfilled. Here, the constraints are obviously regular at any g: Define F : 
L 2 (o,T) ^ K by 

rT 


Then 


Fi9) = [ 9iT)dT-l. 
Jo 

F'{9)h= f h{T)dT, 

Jo 


and hence F'{g) : L^(0,T) —> K is surjective for all g € L^(0,T). 

A simple pointwise discussion of (3^ leads to the following complementar¬ 
ity conditions for almost alH € [ 0 , T : 


0 if i^g(t) + fj- + K / / <p(x, s)u(x, s — t) dxds > 0 
9{t) = { (31) 

/3 if I'git) + n + K J J ip{x, s)u{x, s — t) dxds < 0 . 

If the expression in right-hand side above vanishes, then we obviously have 

g{t) = - ^p + K JJ (p{x, s)u(x, s — t) dxds^ . 

In a known way, the last three relations can be equivalently expressed by the 
projection formula 

9{t) = IP[o./ 3 ] ^ Jl + K JJ^ Cp{s)u[x, s-t) dxds^ ^ , 


where P[o.,fl] : K —>■ [ 0 , /3] is defined by 


P[o.;S](a^) = max(0, min(/3, x)). 
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5 Discussion of (PS) 

Let us now discuss the changes that are needed to establish the necessary 
optimality conditions for the problem (PS) with the particular form (H) of g. 
Now, K,fi, and t 2 are our control variables. Let us denote by U(^K,^ti,t 2 ) the 
unique state associated with (K,ti,t 2 )- 

The existence of the derivatives i = 1,2, and oan 

be shown again by the implicit function theorem. We omit these details, be¬ 
cause one can proceed analogously to the discussion for (PG). To shorten the 
notation, we write 

t = 1,2, 

By implicit differentiation, we find the functions Zi from linearized equations. 
Assume that the derivatives have to be determined at the point (k, ^ 1 ,^ 2 ) and 
fix the associated state u := 'a(K,ti,t 2 ) for a while. Then, zi solves 


{dtZi — Azi + R!{u)zi -\- K zi) (a:, t) 


d_ 

dti 


K 

t2 — ti 



t) dr 


K 

t 2 — ti 


1 

t 2 — ti 



t) dr 


pt 2 

u{x, t — ti) + / Zi{x, t 

Jtl 


t) dr . 
(32) 


Analogously, we find for Z 2 


{dtZ 2 - Az 2 + R'{u)z 2 + KZ 2 ){x,t) = 


— K 

t 2 — ti 


1 

<2 — ti 



t ) dr 


u{x, t — ^ 2 ) 



(33) 


and for Z 3 


{dtZ 3 - Az 3 + R'(u)z 3 + KZ 3 + u)(x, t)= — 


1 


t 2 — ti 

Therefore, the equation for Z 3 is 


nt2 p^2 

/ u{x,t — T)dT + K / Z3{x,t — T) 

Jtl Jtl 


t2-t 
— t) dr 


rt 2 

- / u{x, t — T)dT 
1 Jtl 


K 

{dtZ 3 - Az3 +R'{u)z3 +K.Z 3 ){x,t) - - -— / Z 3 {x,t - t) dr 

C 2 — ti Jt^ 

1 

= - / u{x,t — t) dr — u{x,t). 

h - h J. 


(34) 


Again, we introduce an adjoint equation to set up the optimality conditions. To 
this aim, let (^,^ 1 ,^ 2 ) an arbitrary fixed triplet and M(K,ti,t 2 ) be the associated 
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state function. The adjoint equation is 


{-ipt - Aip + R'+ K(fi) {x,t) 


K 

t2-t 


pt2 

1 Jti 


(p{x, t + t) dr 




dn^p = 0 in if, 


Ud{x,t)) 

in Q, 


(p{x,t) = 0 in 17 X [T, 2T]. 

(35) 

This equation has a unique solution ip € L°°{Q) denoted by <P{K,ti,t 2 ) to in¬ 
dicate the correspondence with (K,ti,t 2 )- Existence and uniqueness can be 
shown in a standard way by the substitution t := T — t that transforms this 
equation to a standard forward equation that can be handled in the same way 
as the state equation. 


Theorem 4 (Derivative of Jg) Let (K,ti,t 2 ) he given, u := he the 

associated state, and ip := ip(K,ti,t 2 ) t/ie associated adjoint state, i.e. the 
unique solution of the adjoint equation ( |35[ ). Write for short 5 := 1/(^2 ~ ti)- 
Then the partial derivatives of Js at {K,ti,t 2 ) are given hy 


dtiJs = >^11 +'^ JJ ip{x,t) ^ u{x,t - t) dr - u{x,t - ti) 

dt2Js = vt2 - jj ip{x,t) j u{x,t - t) dr - u{x,t - t2) 


dxdt. 


dxdt 


dK,Js = vK+ [[ ip(x,t) - [ u{x,t — t) dr — u{x,t) 

J Jq 1° Jti 


dxdt. 


Proof We verify the expression for dtiJs{K,ti,t 2 ), the other formulas can be 
shown analogously. To this aim, let Zi = dt^U(^f^di,t 2 ) be the solution of the 
linearized equation (32l. For convenience, we write z := zi within this proof. 
Following the proof of Lemma we multiply (32) by and integrate over Q. 
We obtain 



{dtz ip + V z ■ Vp + {R'(u) -I- k)z p) dxdt 


[[ p{x,t) - [ u(x,t — t) dr — u(x,t 

J Jq Jti 

+ T // / z(x,t — t) drdxdt. 

0 JJo Ju 


ti) dxdt 


(36) 
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Next, we multiply the adjoint equation (35l by z and integrate over Q to find 


JJ {—zdt(fi + '^z-\/(p+{R'{u) + K)zip)dxdt 


Q ^ti 


(p{x,t + t) dr z{x,t) dxdt + JJ cq{u — Ud) z dxdt. 


(37) 


Now recall that 


1 


0, elsewhere. 


Therefore, we can write 
rt2 


// (p{x,t) / z{x,t — t) drdxdt = // / ip{x,t)g{T)z{x,t — t) drdxdt 

JJq Jti JJqJo 

'll I 9{T)‘Pix,t + T)dTz{x,t)dxdt = - jj j (p{x,t + T)dTz{x,t)dxdt, 

J Jq Jo ° J Jq Jti 


where the second equation follows from (291. In view of 


Q 


{—z) dt<f dxdt = // Lpdtzdxdt, 


Q 


a comparison of the equations (36) and (371 yields 
cq{u — Ud) z dxdt = 


K 

I 


[[ if(x,t) - [ u{x,t — t) dr — u{x,t — ti) 

JJo Ju 


dxdt. 


Since 


^ IIq ~ Udfdxdt = JJ^ CQ{u(^^di,t2) - Ud) Zi dxdt 


the first claim of the theorem follows immediately. 


□ 


As a direct consequence of the theorem on the derivative of Js, we obtain 
the following corollary. 

Corollary 2 (Necessary optimality condition for (PS)) Let {R,ti,t 2 ) 
be optimal for the problem (PS) and let u := and (p := (p(jidi,t 2 ) 

denote the associated state and adjoint state, respectively. Then, with the gra¬ 
dient V Js{ii,ti,t 2 ) defined by Theorem^with (p and u inserted for tp and u, 
respectively, the variational inequality 

VJs(K,ti,t 2 ) • (k - K,tl - tl,t 2 - ^ 2 )^ > 0 y{K,ti,t 2 ) € Cs (38) 


is satisfied. 
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Since the variable k is unrestricted, the associated part of the variational 
inequality amounts to 

lyR — f I (p{x,t)u{x,t) dxdt + ^^ jj (p{x,t)u{x,t — t 2 ) dxdt = 0. 

J Jq ^2 — ti J Jq 

If (^ 1 ,^ 2 ) belongs to the interior of the admissible set C^, then the associated 
components of V J 5 must vanish as well, hence 




vt2 — 


K 

^2 — ti 
K 



1 

pt2 

/ u{x^t — t) dr — u{x^t — ii) 

Jii 

t2 — ii 


1 

pi2 

/ u{x,t — t) dr — u{x,t — 12 ) 
Jt2 

t2 — ii 


dxdt = 0 


dxdt = 0 . 


Remark 3 (Application of the formal Lagrange technique) To find the form 
of VJs and for establishing the necessary optimality conditions, it might be 
easier to apply the following standard technique that is slightly formal but 
leads to the same result: We set up the Lagrangian function C, 

£{u, K, ti,t2, £) = ^ II cq{u - Udfdxdt - II (p{x, t) ■ 


{dtu — Au + R{u) + Ku){x,t) -/ u{x,t — T)dT 

h - ti Ju 


dxdt. 


If (K,ti,t 2 ) is a given triplet, then the adjoint equation for the adjoint state 
is found by 

duC{u, K, ti,t 2 , (p)v = 0 Vu : v{-, t) = 0 for t < 0. 


The derivatives of Js are obtained by associated derivatives of C. For instance, 
we have 

dtiJsiK,ti,t2) = dtiC{u, K,ti,t2,g>) 

if u = inserted after having taken the derivative 

of C with respect to ti. This obviously yields the first two components of V Js 
in Theorem^ For the third component with proceed analogously. 


6 Numerical examples for (PS) 

6.1 Introductory remarks 

The numerical solution of the problems posed above requires techniques that 
are adapted to the desired type of patterns Ud. In this paper, we concentrate 
on numerical examples for the simplified problem (PS), where the kernel g 
is a step function. Although this problem is mathematically equivalent to a 
nonlinear optimization problem in a convex admissible set of the obtained 
patterns are fairly rich and interesting in their own. In a forthcoming paper 
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to be published elsewhere, we will report on the numerical treatment of the 
more general problem (P), where a kernel function g is to be determined. 

In this section, we present some results for the problem (PS), where a 
standard regularized tracking type functional J is to be minimized in the set 
Cs- It will turn out that (PS) is only suitable for tracking desired states Ud 
of simple structure. In all what follows, 17 = (a, 5) is an open interval, i.e. we 
concentrate on the spatially one-dimensional case. 

Compared to many optimal control problems for semilinear parabolic equa¬ 
tions that were investigated in literature, the numerical solution of the prob¬ 
lems posed here is a bit delicate. We are interested in approximating desired 
states Ud that exhibit certain geometrical patterns. If they have a periodic 
structure, then the objective function J may exhibit many local minima with 
very narrow regions of attraction for the convergence of numerical techniques. 
Therefore, the optimization methods should be started in a sufficiently small 
neighborhood around the desired optimal solution. Moreover, the standard 
functional J does not really fit to our needs. We will address the tracking of 
periodic patterns Ud in Section 


6.2 Discretization of the feedback system 


To discretize the feedback equation © , we apply an implicit Euler scheme with 
respect to the time and a finite element approximation by standard piecewise 
linear and continuous ansatz functions (”hat functions”) with respect to the 
space variable. 

For this purpose, the define a time grid by an equidistant partition of [0, T] 
with mesh size r = T/m and node points tj = j t, j = 0,... ,m. Associated 
with this time grid, functions : 17 —>■ K are to be computed that approximate 
u{-,tj), j = 0,.. . ,m, i.e. Uj ^ Based on the functions Uj, we define 

grid functions : Q —>■ R by piecewise linear approximation, 

= -[{tj+i - t)uj{x) + {tj - t)uj+i{x)], if t G [tj,tj+i], j = 0,... ,m. 

T 

By the implicit Euler scheme, the following system of nonlinear equations is 
set up. 


u'^{x, tj+i) — u'^{x, tj) 


- Au^{x,tj+i) -k R{u^{x,tj+i)) 


t2 — ti 


pt2 

/ u'^ {x^tjj^i — s) ds — Ku'^ {x,tj + i), j = 0,.. 

Jti 


m — 1. 


The spatial approximation is based on an equidistant partition of 17 = (a, b) 
with mesh size h > 0. Here, we define standard piecewise affine and continuous 
ansatz functions (hat functions) t/fi : 17 —^ R and approximate the grid function 
u” by ul, 

n 

ul{x,tj) = ^Uj*(/)i(a;) 

2 = 0 
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with unknown coefficients Ujt G K, j = 0, ..., m, t = 0, ..., n. 

To set up the discrete system, we define vectors u-^ G by = 

(ttjo, • • •, j = 0,..., TO. Moreover, we establish the mass and stiffness 

matrices 

M := ( f (j)kix)(j)e{x) dx\ , ^ ■= ( [ 

\Jo J k,e=o \J 

and, for j = 0 ,..., to, the vectors 

i?(u^) ■= ( [ Riul{x,tj)) (j)i{x) dx\ 

V*/ f2 J i—O 

K{u^) ■■= (/^/ {ul{x,tj)-s)ds 

Remark 4 To compute the integrals R{u^), we invoke a 4 point Gauss integra¬ 
tion. Notice that the functions x i—)■ R{ul^{x,tj)) are third order polynomials 
so that the integrand of i?(u^ ) is a polynomial of order 4. Here, the 4 point 
Gauss integration is exact. 

For the computation of the vectors K{u^), we use the trapezoidal rule. 
Here, for tj — s < 0, the values uo{x, tj — s) must be inserted. To increase the 
precision, the primitive of Uq is used. The complete discrete system is 

Mu-^'^^—Mu-^+tAu^~^^+tR(u^~^^) = j = 0,1 ,.. .to—1 . 

We define 

F’(u) := [(1 -I- tk)M + tH]u -|- r i?(u) — tK{u). 

In each time step, we solve the nonlinear equation F'(u) = MW to obtain 
ui+i^ To this aim, we apply a fixed point iteration. 


i:)(j)'^{x) dx 


k,e=o 


i>i{x) dx 


i=0 


6.3 Numerical examples for the standard tracking functional 

For testing our numerical method, we generate the desired state Ud as solution 
of the feedback system for a given triplet {k,ti,t 2 ), he. 

If the regularization parameter v is small, then the numerical solution of (PS) 
should return a vector that is close to {k,ti,t 2 )- 

In all of our computational examples, we selected a small number d > 0 
so that the restriction t 2 — ti > 6 was never active. Moreover, the Tikhonov 
regularization parameter ly was set to zero, because a regularization was not 
needed for having numerical stability. 

In the first numerical examples of this paragraph, the aim is to generate 
desired wave type solutions that expand with a given velocity. 
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Example 1 (Desired traveling wave with pre-computed Ud) We select f2 = 
(—20,20), T = 40, Ml = 0, M 2 = 0.25, U 3 = 1, k = 0.5. Moreover, we take 


as initial function 
uo{x,t) 


1 

2 


— tanh 


X — vt\ 

2V2; 


a; G 17, t < 0 , 


where v = (1 — 2 m 2)/-\/2 is the velocity of the uncontrolled traveling wave given 
by Mo- Following the strategy explained above, we fix the triplet (k,ti,t 2 ) by 
(0.5,0.456,0.541) and obtain Ud = '^(Ktit 2 )- This is a traveling wave with 
a smaller velocity Vd ~ 0.25 due to the control term. To test our method, 
we apply our optimization algorithm to find (k, 11 ,^ 2 ) such that the associated 
state function u coincides with Ud- The method should return a result (k, ti, ^ 2 ) 
that is close to the vector (0.5,0.456,0.541). 

To solve the problem (PS) in Example we applied the Matlab code 
fmincon. For the gradient VJs(K,ti,t 2 )) ^ subroutine was implemented by 
our adjoint calculus. In this way, we were able to use the differentiable mode 
of fmincon. 

During the optimization, we fixed k = 0.5 and considered the minimization 
only with respect to (ti,t 2 )- Taking = 0, ^2 = 1 as initial iterate for the 
optimization, fmincon returned ti = 0.456, ^2 = 0.541 as solution; notice that 
we fixed 1 / = 0. The result is displayed in Fig.[T] 



Example 2 (Stopping a traveling wave ) In contrast to Example]^ here we fix 
the desired pattern Ud that is displayed in Fig.|^ left side. This desired pat¬ 
tern was not pre-computed but geometrically designed, i.e. it is a ’’synthetic” 
pattern that shows a traveling wave stopping at the time t « 16. The other 
data are selected as in Example 

In the optimization process for Example]^ we fixed k = —1.5. The initial 
iterate was ti = 0, ^2 = 1; fmincon returned ti = 0.05, ^2 = 0.94 as solution. 
The optimal state is displayed in Fig.[^ 
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The next example shows that the applicability of the standard tracking 
functional J of (PS) is limited to simple patterns u^, e.g. wave type solutions 
of constant velocity. 


Example 3 (Periodic pattern) Also here, is a synthetic pattern that was 
not precomputed. In f? = (—20, 20), we define 

Ud{x,t) = 3 sin(t - cos(^(a: + 40))). 

Notice that this function Ud obeys the homogeneous Neumann boundary con¬ 
ditions. It is displayed in Fig.[^ left. Since such a periodic pattern cannot be 
expected for small times, we consider the tracking only on [T/2, T]. Therefore, 
here we re-define the objective functional J by 

Js{K,ti,t2) := [ [ {u(^^di,t2)-'^d)'^d.xdt+^{K^+tl+tl). 

Jt/ 2 Jq ^ 


During the optimization run, we fixed the values k = —2 and = 0 and 
optimized only with respect to f 2 - Starting from <2 = 2, the code fmincon 
returned the solution t 2 = 3.71. At this point, the Euclidean norm of VJs is 
|VJs(—2,0,3.71)1 « 0.045. This is a fairly good value and indicates that the 
result should be close to a local minimum. Nevertheless, the computed optimal 
objective value is very large, 

Js(-2,0,3.71) = 1.32 • 10^ 


Remark 5 In this and in the next examples, we fix ti = 0. We observed in 
our computational examples that the optimization with respect to {K,t 2 ) and 
^1 = 0 yields sufficiently good results. Moreover, we found examples, where we 
got the same optimal objective value of J for very different triplets (K,ti,t 2 )- 

The computed optimal state for Example is far from the desired one. 
In particular, the temporal periods are very different. The reason is that the 
standard quadratic tracking type functional J is not able to resolve the desired 
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Fig. 3 


Example 


i 


Desired periodic pattern (left) and computed optimal state u (right). 


periodicity. The main point is that the L^-norm of the difference of a time- 
periodic function 1 1 —>■ Udit) and its phase shifted function t ^ Udit — s) can be 
very large, although both functions have the same time period. For instance, 
in Q = (—20, 20) x (0,40) we have 



-cos(^(a;-h20))) 


3sin(t — 3 


cos(^(a;-h20)))) 


2 

dxdt = I, 


where I ~ 1.4374-10^. This brought us to considering another type of objective 
functionals that is discussed in the next section. 


7 Minimizing a cross correlation type functional 

The cross correlation 


As Example showed, we need an objective functional that better expresses 
our goal of approximating periodic structures. This is the cross correlation 
between u and Ud- Moreover, in the functional, we have to observe a later part 
of the time interval, where u already has developed its periodicity. 

Let us briefly explain the meaning of the cross correlation. Assume that 
iCd : R —t K is a periodic function and x ■. [0, T] — >■ K is another function; think 
of a function x that is identical with Xd after a time shift. 

To see, if Xd and x are time shifts of each other, we consider the extremal 
problem 


min 

sGR 


[ {x{t) - Xd{t - s jfdt. 

Jo 


(39) 


If X and Xd are just shifted, then the minimal value in (39 1 should be zero by 
taking the correct time shift s. The functional (391 can be simplified. To this 
aim, we expand the integrand. 


{x{t) — Xd{t — s)Ydt= / x^{t)dt—2 / x{t)xd{t — s)dt-\- / x‘^{t — s) dt. 
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The first integral does not depend on s. Since Xd is a periodic function, also 
the last integral is independent of the shift s. Therefore, instead of minimizing 
the quadratic functional above, we can solve the following problem: 


max 

sGR 


x(t)xd(t - s)dt 
\\x\\L^(0,T)\\xd\\L'^{0,T) ’ 


(40) 


where we additionally introduced a normalization. The result is the so-called 
cross correlation between x and Xd- The largest possible value in (401 is 1; in 
this case, both functions are collinear. 

In the application to our control problems, this induces two equivalent 
objective functionals. The minimization problem (39 \ leads to the optimization 
problem 

( min / {u{x, t) — Ud(x, t — s))^dxdt + ^ +^ 1 +^ 2 )) ■ (41) 

\ Jq 2 / 


mm 

{K,tl,t2)eCs 

The other way around is an equivalent problem that uses the cross corre¬ 
lation, 

(42) 


min Jcorr(K,<l,< 2 ) 

(Kdl,t2)C:Cs 


where 


JcorriK,ti,t2) ■= 1 — max 
sGR 


JJ u{x,t)ud{x,t — s)dxdt 
lkllL2(Q)||Ud||L2(Q) 


+ ^i^^ + tl + 4)- (43) 


For solving (421, we applied the Matlab code pattern search that turned 
out to be fairly efficient in finding global minima for functions that exhibit 
multiple local minima. In the periodic case, we have to deal with multiple lo¬ 
cal minima indeed. Here, the cross correlation based functional Jcorr is more 
useful, as the next figure shows. 


Example 4 (Multiple local minima) We pre-compute the desired state by Ud = 
■“(-2,0,2.5) and consider the functions k i - a - Js{k, 0, 2.5) and k g ->• Jcorr{K, 0, 2.5) 
for K around the (optimal) parameter R = —2. 


The two functions defined in Example]^ are shown in Fig.|^ The function 
K G->■ J 5 («:, 0,2.5) behaves more wildly around k = —2 than the function 
K GA- Jcorrin, 0, 2.5) that is based upon the cross correlation. 

Let us re-consider the optimization problem of Example but now by 
the cross correlation based optimization problem (42l. Here, we apply the 
following strategy: We keep = 0 fixed and optimize only with respect to 
(k, ^ 2 )- Moreover, at the beginning we fixed n = —2 and optimized with respect 
to K in a second run. The computed solution was ^2 = 2.7631 with a value 
Jcorr = 0.1604, (ti = 0, K = —2); as before the Tikhonov parameter was 
selected as z/ = 0. The computed u of this first step is shown in Fig.[^ Now 
the agreement, in particular of the temporal period, is much better. 









28 


P. Nestler et al. 




Fig. 4 Example^ goal functions k i— > Jg(K, 0,2.5) (left) and k i— > Jcorr(K, 0, 2.5) (right). 


Next, we performed an alternating search for (k, ^ 2 ) starting with the result 
obtained in the first step. We obtained t 2 = 2.7631, k = —2.4318 and the 
improved objective value Jcorr = 0.1229. This improvement is graphically 
hardly to distinct from Fig. 



Fig. 5 Example 1^ Desired pattern (left) and computed pattern after minimizing Jcorr 
with respect to t 2 (right) for fixed k = —2, ti = 0. 


Finally, we consider another example with synthetic Ud that has a larger 
period than Ud of Example 

Example 5 We consider again fl = (—20, 20) and observe Ud only in the time 
interval [20,40]. For Ud we select 


Ud{x, t) = 3 sin 



Again, ti = 0 is fixed and the iteration is started with ^2 = 2, k = —2. 
The optimal control parameters are ^2 = 6.94191, k = —2.28 with computed 
optimal objective value Jcorr = 0.1209. The computed optimal state is shown 
in Fig.[^ 
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